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Abstract. Photoassociation of ultracold atoms induced by chirped picosecond pulses is analyzed in a non- 
perturbative treatment by following the wavepackets dynamics on the ground and excited surfaces. The 
initial state is described by a Boltzmann distribution of continuum scattering states. The chosen example 
is photoassociation of cesium atoms at temperature T=54 fiK from the a'^E^{6s, 6s) continuum to bound 
levels in the external well of the 0^ (6s + 6^3/2) potential. We study how the modification of the pulse 
characteristics (carrier frequency, duration, linear chirp rate and intensity) can enhance the number of 
photoassociated molecules and suggest ways of optimizing the production of stable molecules. 

PACS. 33.80.Ps Optical cooling of molecules; trapping; 33.80 -b Photon interactions with molecules 33.90 
+h New topics in molecular properties, interaction with photons 33.80.Gj diffuse spectra; predissociation, 
photodissociation 

1 Introduction and to create excited molecules via the photoassociation 

reaction [S]; subsequently, these molecules are stabilized, 
The various routes leading to the formation of cold and by spontaneous emission or other radiative coupling, into 
ultracold molecules are actively explored Hj. Non-optical bound vibrational levels of the ground electronic state 
techniques like buffer gas cooling of molecules and |7|«|iJ|lU|ri|l^]. The translational temperatures thus reached 
Stark deceleration of polar molecules EEI reach temper- are much lower (T < 20 fiK). Nevertheless, such stable 
atures well below 1 K. Another route relies on optical molecules are produced in a superposition of vibrational 
techniques, laser fields being used to cool alkali atoms levels, the most of them very excited. Then, bringing the 



E. Luc-Koenig, M. Vatasescu, and F. Masnou-Seeuws; Cold atoms photoassociation using chirped laser pulses 



t=t. 



molecules to the lower vibrational level (v=0) , thus mak- 
ing vibrationally cold molecules, is an important issue. 

Up to now, most photoassociation experiments are us- 
ing continuous lasers, but there are a few papers treating 
the photoassociation with pulsed lasers |13I14I15| . Theo- 
retically, several time-dependent studies of photoassocia- 
tion in the ultracold regime have been proposed |16I17I18I19I20 | 

Our aim is to investigate the possibility to control cold 
molecules formation by use of chirped laser pulses. In- 
deed, it was shown that picosecond frequency-swept laser 
pulses produce more selective excitation and better pop- 
ulation transfer than transform- limited pulses with the Fig. 1. Scheme of the photoassociation process with a negative 




6s+6s 



same bandwidth, duo to the mechanism of population in- 
version by adiabatic sweeping |21I22I23I24| . 



chirped pulse considered in the present work, illustrated in the 
case of Cs2. The potentials curves correspond to the ground 



continuum state is neglected in the definition of the resonance 
condition. The double well behaviour in the excited curve is a 
particular feature of the chosen symmetry. 



T-,1 , ■ r 1 1 i -ii 1 • 11 1 triplet state a^X'.tfGs -I- 6s) and to 0„ (6s + Qp-i/2) excited 

Photoassociation oi cold atoms with a chirped laser pulse ^ ' ^ ^ f->/^' 

—— electronic state. In the present work, the energy of the initial 
was nrst explored theoretically by Vala et al. [18j, using 

a gaussian packet centered at large interatomic distance 
(i? sa 200 ao) as initial state in the collision of two ce- 
sium atoms at the temperature T = 200 iiK, and show- 
ing that a picosecond pulse can achieve a total transfer representation of the initial continuum state is a station- 
of population under adiabatic following conditions pro- ary coUisional state. Moreover, at a given collision energy, 
posed by Cao, Bardeen, and Wilson |25l26) . Nevertheless, the behaviour of a continuum state varies a lot with the 
such an approach cannot work at collision energies close interatomic distance, and then the distance range where 
to threshold, and if we are interested in the dynamics at the excitation takes place (for example, Rl in Fig. ^| is 
smaller distances (which it is the case if the aim is to pro- of great importance for the photoassociation yield. There- 
duce cold molecules in low vibrational levels), because a fore, a theoretical treatment of the dynamics at very low 
gaussian wavepacket does not adress the actual shape of temperatures needs a quite precise representation of the 
the initial continuum state. Indeed, in an ultracold colli- initial state. 



sion of two atoms, the kinetic energy ksT becomes eas- 



In a previous theoretical paper 12 Oj . we have consid- 



ily smaller than the interaction potential, and the correct ered the photoassociation of two cold cesium atoms with 
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a picosecond chirped pulse that excites several vibrational 
levels in the 0^ (6s + 6^3/2) external well. This reaction 
has been widely studied experimentally with continuous 
lasers, so that the relevant potential curves are well known. 
The time-dependent Schrodinger equation corresponding 
to the two-electronic coupled channels (a'^Z'+(6s-t-6s) and 
0~(6s -|- 6^3/2)) was solved using a propagation method 
|27I28| . and, due to a Mapped Sine Grid representation 
PU] , it was possible to take as the initial state in the pho- 
toassociation a stationary collision state corresponding to 
a temperature of about 50 ^K. The chirped pulse has been 
designed in view of creating a vibrational wavepacket in 
the excited state, which, after the pulse, moves to shorter 
internuclear distances and is focussing at the barrier of 
the double well potential, thus preparing a good initial 
state for the stabilization process into low vibrational lev- 
els of the a?S^{&s + 6s) potential. Besides its property 
of focalization, the specificity of the picosecond chirped 
pulse analyzed in Ref. [201 is that it leads to a strong 
transfer of population inside a spatially limited "photoas- 
sociation window" , and no tranfer outside. Finally, the 
population transfer to the last vibrational levels of the 
ground a^E^{Qs + 6s) state is significant, making stable 
molecules. We have shown that these results can be in- 
terpreted in the framework of a two-state model and in 
the impulsive limit |3L)| . as an adiabatic population inver- 
sion taking place in the "resonance window" swept by the 
pulse. 

In the present work, we consider the same photoassoci- 
ation reaction, and analyse how the results are modified 
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by acting with various pulses, in order to determine how 
characteristics like duration, chirp rate and energy influ- 
ence the photoassociation yield. The results obtained us- 
ing pulses of different durations (ps, tens and hundreds of 
picoseconds) and energies, can then be used to explore var- 
ious regimes of population transfer, nonadiabatic effects, 
and the eventual limits of the impulsive approximation. 
Moreover, as our calculation are performed using a very 
large spatial grid (thousands of oq), the results account for 
the threshold effects intervening in this continuum-bound 
transfer of population taking place at very low energies 
and implying large ranges of distance in the transition (as 
the initial state is a delocalized coUisional continuum). 
The exploratory results thus obtained will be used in the 
calculation of the effects of a second pulse, which can be 
either one identical with the first (to estimate the effect 
of the repetition rate), either a different one (chosen for 
the stabilization of the packet focalized in the 0^^ channel 
to vibrational levels of the a^S^ potential). Obviously, 
this non-perturbative model allows to treat laser pulses of 
strong intensities and non-resonant effects. 

The photoassociation reaction studied here is between 
two cold cesium atoms colliding in the ground state poten- 
tial a'^Z'+(6s -I- 6s), at a temperature T = 54 iiK^ which 
are excited by a pulse laser to form a molecule in a super- 
position of vibrational levels {w} of the electronic potential 
Og (6s + 6^3/2). For a rotational quantum number J=0, it 
can be represented as: 

Cs{<6s^Si/2) + Cs{<6s^Si/2) + nuj{t) 
^ Cs2(0;(6s25i/2 + 6p2p3/2); {v}, J = 0), (1) 
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We consider the excitation by a chirped laser pulse of gaus- 
sian envelope, having a time-dependent frequency cj(t)/27r 
which varies linearly around the central frequency LOi^jlix 
reached at t = tp, and red-detuned by 5£* relative to the 
Z?2 atomic resonance line: 

huj(ip) = hujL = ?iwat - (2) 

hiOat is the energy of the atomic transition 6s 6p3/2. 
The central frequency of the pulse, lj L/2Tr, or equivalently 
the detuning determines the crossing point Rl of the 
two electronic potentials dressed by the photon with the 
energy Tilol (see Fig. In the present calculations the 
detuning is fixed at (52*=2.656 cm^^, corresponding to 
Rl = 93.7 ao and to resonant excitation at t = tp of 
the level wo=98 in the external well of the 0^ (6s + 6^3/2) 
potential. Then, taking the origin of energy as being the 
dissociation limit 6s -|- 6^3/2 of the 0^ (6s 4- 6^3/2) poten- 
tial, the binding energy of the wq level is E'u,, = In 
this paper the origin of energy is chosen at the dissocia- 
tion limit (6s + 6s) of the a^S^ potential, the energy of 
the Og vq level being then equal to 0. 
The potentials used in our calculation have been described 
in a precedent paper |2Q|. The outer well of the 0^ (6s -|- 
6^3/2) excited potential was fitted to photoassociation spec- 
tra by Amiot et al \6\\ and matched to ah initio calcula- 
tions at short and intermediate range ■ The a'^Z'+(6s, 6s) 
potential has been chosen in order to reproduce correctly 
the scattering length L w 525 oq and the asymptotic be- 
haviour —Cq/R^ with Cq— 6828 au [3T] (the short range 
part extracted from the Ref. ^33: being slightly modified 
for that purpose). 
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The paper is organized as follows: in Section 2 we 
discuss the representation of the initial state in the pho- 
toassociation dynamics. Section 3 describes the numerical 
methods used to perform time-dependent propagation on 
very large spatial grids. Section 4 shows the two-channel 
model used for the photoassociation, the basic formulae for 
gaussian pulses with linear chirp, the choice of the pulse 
parameters, and treats specific features of the excitation 
with a chirped pulse, as the energy range excited reso- 
nantly or not resonantly, and the condition for an "adi- 
abaticity window" during the pulse duration. Section 5 
compares the results obtained with pulses differing as du- 
ration, energy, and chirp. In Section 6 we show the calcu- 
lation of the photoassociation probability from a thermal 
average over the incident kinetic energies, and of the num- 
ber of molecules photoassociated per pump pulse. Section 
7 presents a discussion on the possible ways for the opti- 
mization of the process. Section 8 is the conclusion. The 
article has three appendix, containing: Appendix A (Sec- 
tion 9) - the energy normalization of the ground state 
continuum wavefunction calculated in a box; Appendix B 
(Section 10) - the calculation of the chirp rate in the time 
domain for focussing the excited vibrational wavepacket 
at the outer turning point; Appendix C (Section 11) - the 
overlap of the initial continuum with the 0~ vibrational 
wavefunctions. 
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2 Representation of the initial state in the 
photoassociation dynamics 

We shall briefly review here the description of the initial 
state of the photoassociation process, as a thermal equi- 
librium state in a gas of cold alkali atoms, lying in a trap 
of infinite dimension with the laser turned off (then not 
interacting with the electromagnetic field), and being ini- 
tially in their {ns) ground state. Supposing low atomic 
densities, which allows to consider the binary events iso- 
lated enough, the N-body {N > 3) interactions can be 
ignored, and the Hamiltonian will be reduced to a sum- 
mation of non-interacting pair Hamiltonians, treated in 
the Born-Oppcnheimer approximation 34 . Then, for a 
colliding pair of atoms, the molecular Hamiltonian 
of the relative motion can be separated from the center 
of mass motion which does not play any role in photoas- 
sociation. Assuming thermal equilibrium at temperature 
T, the initial state is a statistical mixture of eigenstates 
corresponding to Hm, being described by the density op- 
erator: 

P=^^, P=l/{kBT) (3) 

with Z = Tr{e^^^'^'} the partition function (fc^ being the 
Boltzmann constant). For a gas of low density composed of 
non-interacting pairs of atoms in a volume V, the partition 
function is |35| : 

z.omv, o,T,= <™!^!, (4, 

/i being the reduced mass of the diatom. 

The density operator ijjjl can be decomposed into par- 
tial waves I describing the relative rotational angular mo- 
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menta: 

p=lY.i2l + l)pl pi = e-P"" (5) 

where the molecular Hamiltonian Hi ~ ^ + ^g{R) + 
^^-jj^j^- R designs the interatomic coordinate, = T 
is the kinetic energy operator, and Vg{R) is the potential 
of the ground electronic state (neglecting the hyperfine 
structure, g represents one of the states X^S'^{ns + ns) 
or a^S:^{ns + ns)). We shall restrict our study a,t I — 
(s wave), which is a good approximation for low energy 
collisions, and show here calculations for cesium photoas- 
sociation, taking a'^Z'+(6s-f 6s) as electronic ground state. 
Taking the dissociation limit (6s + 6s) of the a'^Z'+ poten- 
tial as origin of energy, the initial state of the photoas- 
sociation process will be described by a density operator 

Ps = ^j^ dEe-^^'lE X E\ (6) 

where \E >= P^^, I — 0, E > are continuum eigenstates 
(normalized per unit energy) of the Hamiltonian Hi^q, of 
energy E and corresponding to the relative motion in the 
electronic potential and to Z = (s -wave). 

In the following wc shall discuss the representation of 
the initial state in the study of the photoassociation dy- 
namics by time-dependent propagation. An evaluation of 
the initial density matrix by gaussian wavepackets can be 
valid for sufficiently low laser detunings (very large dis- 
tance R) and sufficiently high temperatures T jl8ll9| . On 
the contrary, at very low temperatures and for the pho- 
toassociation in more profound vibrational states of the 
excited molecular potential, the potential energy Vg{R) 
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becomes significant relative to the kinetic energy term, 
and the calculation needs to account for the nodal struc- 
ture of the ground state. In this case, the evaluation of 
the initial density matrix by stationary collision states is 
required. In this paper we present calculations correspond- 
ing to this second case, showing that such delocalized con- 
tinuum states can be represented numerically and used in 
time dependent propagation. 

2.1 Evaluation of the initial density matrix by gaussian 
wavepackets 

For a sufficiently small detuning S"^, the crossing point 
Rl of the electronic potentials dressed by Tilul occurs at 
large distance. In the case of the excitation with a contin- 
uous laser, and for a wavepacket localized around R^, the 
crossing point Rl completely determines the R-domain 
governing the photoassociation process. For the photoas- 
sociation with a pulse, other factors strongly influence the 
distance range which can be excited: the duration and the 
intensity of the pulse, its chirp rate (see Sec. 14.4(1 . and the 
dynamics depends on the spatial localization or delocal- 
ization of the initial state. 

For sufficiently high temperatures T, one can expect 
that the kinetic energy dominates the interaction potential 
Vg{R) ~ — at large R distances: 

kBT:^^, (7) 

resulting in a relatively uniform R-dependence for the 
density matrix. If the inequality iQ is valid, the den- 
sity operator can be decomposed in a set of incoming 
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gaussian wavepackets with mean velocity vq equal to the 
most probable velocity of a Maxwell-Boltzmann distri- 
bution (then the corresponding momentum po = fivg = 
\/2fikBT — V < Po > 5 ^-iid the momentum spread is de- 
fined as Apfi ~ The width AR of the radial density 
distribution is related to the momentum spread for the 
relative motion by ARAp^ w |, then depending on the 
temperature jl8ll9| . 

In the case of the a?'S^{Qs + 6s) state of the Cs2, for 
the temperature T « 50 pK studied in the present paper, 
the relation lO is valid for R ^ 185 ao, which corresponds 
to resonant excitation of vibrational states with v 3> 130 
in the 0^ (6s + 6733/2) external well, or to a laser detuning 
(^2* < 0.35 cm~^ At T « 50 fiK (a collision with rela- 
tive velocity vq « 10 cm/s), the width corresponding to 
a radial wavepacket is very large: AR « 160 ao (in cold 
collisions the large de Broglie wavelength describing the 
relative motion is the sign of the wavefunction delocaliza- 
tion). 

Gaussian wavepackets have been succesfully used in 
the study of photoassociation [TS| at a large distance Rl = 
200 ao and a temperature T = 200 iiK, for a pulse dura- 
tion (75 ps) sufficiently small to avoid observable motion 
on the ground surface. 

On the opposite, in a precedent paper jj^j studying 
the photoassociation with a continuous laser in the elec- 
tronic state lg(6s -I- 6^3/2) of Cs2, at Rl — 90 ao and for 
a temperature T — 125 p,K, we have shown that a gaus- 
sian wavepacket having a large radial width and a small 
momentum pq is spreading more rapidly than it is mov- 
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ing; being quite difficult to gain insight on the relevant form, but if one has to account for a range of distances 

dynamics by using it, a possible choice is to consider a such as Rl ^ Rn it is necessary to describe the initial 

much smaller spatial width, implying a bigger momentum density matrix as an incoherent average over a thermal 

width. Components with much bigger momenta are then distribution of stationary energy-normalized continuum 

introduced in the representation of the initial state, falsi- wavefunctions \E >, as in Eq. ©. Each incident compo- 

fying the dynamics in the corresponding channel. nent will contribute incoherently to the density probabil- 

Therefore the representation of the initial state by gaus- ity I'Z'q- (R, t)\'^ transferred at the time t in the 0^ surface, 

sian wavepackets cannot be used at very low temperatures with a weight equal to ^—g— . 
if we are interested in the real dynamics at small distances 

and in the initial molecular channel, because it hardly 2.3 Representation of the energy-normalized initial 

adresses properly the physics at distances different from continuum state on a spatial grid. 

Rl. 

In our method, the eigenstates of the Hamiltonian asso- 
ciated with the ground electronic channel are calculated 

2.2 Evaluation of the initial density matrix by 

through the Sine Grid Representation |29l28j . in a large 

stationary collisional eigenstates 

box of radius Lfj which introduces a discretization of the 

At very low temperatures, the kinetic energy fc^T can continuum (only continuum states having a node at the 

easily become smaller than the potential energy ^, even boundary of the box are obtained) and supplies wavefunc- 

for large values of R. Then it is absolutely necessary to tions normalized to unity in the box. We shall label > 

account for the nodal structure of the ground state wave- the energies of these states and (i?) the corresponding 

function, which is responsible for the intensity minima ob- wavefunctions having n nodes in the range Rn < R < Lr 

served in the photoassociation spectra. This is the case for > integer). The energy resolution in a box of length 

the studied example of photoassociation in the 0^ (6s + is determined by the energy difference SE correspond- 

6P3/2) channel in Csa, at a detuning 61* « 2.4 cm"\ ing to neighbouring eigenstates of the box: 



i?i w 94 an, and T w 50 uK. Indeed, in the ground r^j-^ t /„ ^ t2 2 

dE Tin 2E„ n n 

OE \e^^ \e=e„^ -^W = —jT'^ (8) 

state a^Z'+ (6s + 6s), all the stationary wavefunctions in- ^B. \j M f^^R 

volved in the presently studied energy range have the same Then a large box is necessary to have a sufficiently high en- 
nodal structure until Rn — 82.3 ao (the position of the ergy resolution (small 5E), able to represent low energies 
last common node), differing only by a normalization fac- corresponding to the studied temperature T « 50 iJ,K. For 
tor (see Fig. |2K) )• For photoassociation at large distance the present calculations we use a box with the width = 
Rl S> Rn the distribution of the contimmm states is uni- 19250 ao, in which the continuum state having the energy 
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Fig. 2. a) Stationary continuum states in the o^X'+(6s + 6s) 
potential, corresponding to energies E = ksT, which are cal- 
culated on a grid of length Lr — 19250 ao and normalized 
at 1 in the box. Rn is the last node common to the coUisional 
wavefunctions in the studied domain of energies. With the con- 
tinuum thick line is represented the continuum wavefunction 
presently chosen as initial state in photoassocciation at the 
temperature Eo/ks — 54 fiK. Other continuum wavefunctions 
obtained in the box are represented, corresponding to the fol- 
lowing values of E/kB- 36.6 nK (dot-dashed line), 0.33 ^iK 
(short-dashed line), 0.11 mK (long-dashed line), and 0.32 mK 
(thin continuum line). The last bound state v" = 53 of the 
triplet potential is also shown (with the wavefunction dimin- 
ished by a factor of 5) ; b) The wavefunction corresponding to 
the level vo = 98 in the outer well of the Og (6s -f- 6^3/2) poten- 
tial, which is excited aX t = tp {Ey^ = —2.656 cm~^ under the 
dissociation limit 6s -|- 6^3/2) • 



En = ksT, with T = 54 iiK, is the 38th state in the cor- 
responding discretized continuum. This leads to a resolu- 
tion at threshold |^ = 30 nK and |E/feB=54/iX = 2.6 
^.K for the continuum description around 54 iiK. Such a 
large box allows to account correctly for the threshold be- 
haviour in the presently studied cold collision process and 
to perform explicitely an average on the thermal energy 
distribution. 

For an energy resolution 5E, the longest relevant time 
scale which can be studied in the problem is Tmax = h/5E. 
Any energy spacing smaller than 5E is not resolvable in 
the box of length Lr (for example, an eigenstate charac- 
terized by a vibrational period bigger than Tmax cannot be 
distinguished from the continuum positive energy states in 
the box determined by this Ln) . In the present case Tmax 
is 2.9 /is, being larger than the radiative lifetime in the 0~ 
state, Trad ~ 15 ns; but the spontaneous emission is not 
introduced in the present problem, where the dynamics is 
followed for times smaller than Trad- 

The representation of the initial density matrix by sta- 
tionnary continuum states needs a normalization per unit 
energy of the continuum wavefunctions. The basic formu- 
lae are shown in the Appendix A: one shows that the re- 
lation between the wavefunctions normalized per 
unit energy and the wavefunctions (j)E„{R) normalized to 
unity in the box is the following: 



f) E 



(9) 



In the presently Grid representation the density of states 
is calculated explicitely from the energy spacing between 
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two neighbouring states of the discretized continuum: 3.1 Spatial representation of the wavefunctions 
dE 

\e=e = En+i — En (10) The radial dependence of the wavepackets g{R, t) prop- 

dn 

agating on both surfaces is represented using Mapped Grid 
Methods j36l29j . The mapping works with a change of 

3 Numerical methods 

variable (from R to the adaptive coordinate x) made as 
The dynamics of the photoassociation process is performed to take into account the variation of the local de Broglie 
by solving numerically the time-dependent Schrodinger wavelength A{E,R) (see Eq. as a function of the in- 
equation for the ground Vg (R) and excited Ve (i?) poten- ternuclear distance R. Namely, the Jacobian of the trans- 
tials coupled by the electromagnetic field, and propagat- formation is chosen proportional to the local dc Broglie 
ing the wavepackets on the ground and excited surfaces, wavelength: J{x) — dR/dx — f3A{Emax, R), where E^ax 
The cold collision of the two atoms on the ground surface is the maximum energy involved in the problem. This 
Vg (i?) is generally characterized by a huge local de Broglie leads to the implementation of a large spatial grid by 
wavelength yl(i?, i?), as the relative momentum pfl is very using only a small number of points. For example, in 
low: our calculations we use N=1023 points for a spatial ex- 

7 2 tension ~ 19250 ao, describing uniformely regions of 

A{E, R) = ^ (11) 

■\/2/i[ii Vg{R)\ PR space where A{Emaxi R) varies by several orders of magni- 

As we have shown in the precedent section, we choose tude. The use of an enveloping potential yem)(^) (equal or 

as initial condition a stationary continuum state of the deeper than the crossing Ve and Vg potentials dressed by 

ground surface of very low energy corresponding to a tem- the field) allows the definition of a common Jacobian and 

perature T = 54 ^K. A grid of very large extension is therefore a single a;-grid to describe both surfaces. In the 

needed to represent correctly such an initial state and new coordinate a;, the kinetic energy operator expresses 

the last bound states of the ground potential which are simply as product of operators J(a;)~^/^ and d/dx, with- 

populated during the photoassociation process, and whose out introducing derivatives of J{x) i36l29| . The potential 

wavefunctions extend at large interatomic distances. We energy operator reduces simply to V{x). 



consider a grid extending from Lq (a distance slightly A collocation method |36I29| is used to define the rep- 
smaller than the repulsive walls in the potentials, here resentation of the wavefunctions in the N points Xi of 
Lq—S slq ) to Lii = 19250 ao in the present paper. This the grid (x). This supposes that any wavefunction ^(x) 
very large grid allows to represent correctly the wavefunc- is expanded on a set of N basis functions. Instead of the 
tions implied in the process and to account for the thresh- usual plane wave expansion (Fourier expansion), we use 
old behaviours at low energies in the realistic potentials. the Sine expansion recently introduced by Willner et al. 
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AH basis sine functions have nodes at the bound- 
aries Lq and Lfl of the grid, which, by choosing a suffi- 
ciently small /3- value 129] , permits the suppression of the 
so-called "ghost" levels appearing in the solution of the 
stationary Schrodinger equation and then susceptible to 
falsify the dynamics of the system. We have verified that 
a value (3 = 0.52 is necessary for avoiding the appearance 
of "ghost" levels. 

To determine the stationary eigenstates of the ground 
or excited potential (which are used either for the selection 
of the initial state, or to analyze the wavepackets evolution 
in terms of their decompositions on stationary states) an 
auxiliary cosine basis set is introduced, allowing to eval- 
uate analytically the first order derivatives d/dx involved 
in the kinetic operator |29| . 

In time propagation calculations one has to multiply 
repeatedly the initial state wavefunction by the Hamilto- 
nian matrix. Therefore one has to express the kinetic en- 
ergy operator in terms of discrete sine and cosine Fourier 
transformations 37 , for which we use the efficient fast 
Fourier algorithms. This imposes = 2^ — 1 {p integer) 
for the number of grid points (N=:1023 for the present 
calculation) . 

Lastly, in time propagation calculation the initial state 
is normalized to unit in the box of size Lji. Therefore the 
populations -Pg- (0 or P3^^{t) calculated at a time t are 
depending on the value chosen for L^. Nevertheless, due to 
the normalization condition P^-{t) + P3^^^[t) = 1, which 
is verified for every t, the calculations give directly the 
relative population on each surface. 
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3.2 Time Evolution 

The time-dependent Schrodinger equation for the two cou- 
pled channels (see Eq. (llTfl ) is solved by propagating the 
initial wavefunction using a Chebychev expansion of the 
evolution operator exp[—iflt/h] |27I28| . Propagation is 
performed in discrete steps with a time increase At much 
shorter than the characteristic times of the problem (pulse 
duration, vibrational periods, Rabi periods). In the present 
study of photoassociation with chirped pulses the Hamil- 
tonian H{t) is explicitely time-dependent, and a discrete 
description of H{t) is introduced. During the step corre- 
sponding to time propagation from ti to ^2 = ^1 + At, 
the Hamiltonian is supposed to be time-independent and 
chosen as H{tm), where t„i = {ti + t2)/2. Such a proce- 
dure introduces errors of the order of magnitude (At)^. 
Since presently the initial state is a stationnary collisional 
eigenstate in the potential, instead of a gaussian 

wavepacket, the dynamics results only from coupling with 
the laser pulse, but not from the evolution motion charac- 
teristic for a nonstationary state. Moreover, the dynamics 
is studied within a very large box. Consequently, even af- 
ter a very long propagation duration of 15 ns (the order 
of magnitude of the spontaneous emission time for the vi- 
brational levels excited in the 0^ state), the wavepackets 
dynamics in the range of distances relevant for our prob- 
lem (in which the hyperfine structure is neglected) is not 
influenced by the external boundary of the box. There- 
fore there is no reflection of the wavepackets at L/j and 
then it is not necessary to define outgoing wave boundary 
conditions for the wavepackets (either by transferring the 
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outgoing part to another grid or by introducing an 
imaginary absorbing potential [5^). 

In the present problem, using a grid with N=1023 
points, which for /3 ~ 0.52 gives a box of extension Lj^ = 
19250 ao, the energy resolution at the '^S^ dissociation 
threshold is determined by SE/ks ~ 100 nK. For a grid 
with only 511 points, the extension is reduced at Lji = 
2450 ao (about 7.8 times smaller than Lu) and the en- 
ergy resolution is strongly diminished, as SE/ks ~ 6 fiK. 
Since we study the photoassociation taking an initial sta- 
tionary state corresponding to T~ 50 fiK, a grid of large 
extension is needed to attain the necessary resolution in 
the representation of the threshold processes. As we have 
shown before, this is possible due to the mapping proce- 
dure, which drastically reduces the number of grid points, 
and due to the implementation of the Sine Basis represen- 
tation, which eliminates the participation of the "ghost 
levels" in the dynamics. 

For the step At « 0.05 ps used in our time propaga- 
tion, H!P' is calculated 112 times during each interval At. 
Therefore the implementation of the fast Fourier trans- 
formations to calculate T'F avoids prohibitive calculation 
times, allowing in principle to study the evolution of the 
wavepackets during time durations as long as 10 ns (which 
corresponds to a repetition rate of 10^ Hz for the photoas- 
sociating laser pulse) . Obviously, such long evolution times 
need to introduce in the problem other physical processes, 
as the spontaneous emission {Trad = 15 ns), which is not 
taken into account for the shorter times of our present 
calculations. 
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4 Two-channel model for the 
photoassociation with a gaussian chirped 
pulse 

4.1 Basic formulae for gaussian pulses with linear chirp 

We shall briefly review the main relations concerning pulses 
with gaussian envelope and linear chirp [25126120] . The 

electric field has the amplitude So, a gaussian envelope 
f{t), the carrier frequency 10^/211^ and a phase ^p{t): 

£{t)^£af{t)cos[ujLt + v{t% (12) 

The gaussian envelope f{t) is maximum sXt — tp and has 
the full width at half maximum (FWHM) equal to V2tc. 

/(t) = ./^exph21n2(^)2], (13) 

V Tc Tc 

The duration tl characterizes the spectral width 6lu~ 
AIh2/tl in the frequency domain: indeed, £{w), which is 
the Fourier transform of £{t), displays a gaussian profile 
with FWHM equal to V^Slo. 

The ratio ^ > 1 characterizes the chirp: the value 
Tc = tl corresponds to an unchirped transform-limited 
pulse with duration r^, which is streched to tc by chirp- 
ing. The linear chirp is determined by the chirp rates: x 
in the time domain (calculated as the second derivative of 
(p{t), the phase of the field £{t)), and in the frequency 
domain (the second derivative of the phase of £{lo)): 

The choice of and x imposes the pulse duration tq (see 
Eq. 114|l l. and then the chirp rate <l>" in the frequency 
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domain: 



'C'L 



(15) 



(4 In 2)2 

Among the four parameters tc, ti^, x and <P" , only two 
are independent (see Eas. 1141 and 1151) . Another indepen- 
dent parameter is the coupling Wl- Then, the choice of 
three independent parameters will determine the nature of 
the transfer between the two channels, adiabatic or non- 
adiabatic, as it will be shown in Sec. 14.51 

Chirping a pulse increases its duration and decreases 
its maximum amplitude: £m ^ ^o^^J^ 1^ £o- As a result, 
the chirp does not change the energy Epuise carried by the 
field, which is proportional to the square of the amplitude 
£q and to the temporal width of the transform limited 
pulse: Epuise^^ j^^\£{t)\''dt^ ^V^' ^it^ = 

ceo c2 
2 ^0- 

The instantaneous frequency Lo{t)/2'iT is given by the 
derivative of the rapidly oscillating term in £{t): 



Lj{t)^LJL + ^^UJL+x{t~tp), (16) 

at 



and varies linearly around the central frequency u^i^jlix. 

4.2 The two-channel coupled equations 

We study the dynamics in the ground and excited states 
coupled by the pulse described in Ea. H12(l . for which we 
consider linear polarization Bl. In the dipole approxima- 
tion, the coupling term between the two electronic chan- 
nels is written as: — Dge(_R) • e]_,£(t) w Dg^£{t), where 
Dge(-R) is the R-dependent matrix element of the dipole 
moment operator between the ground and the excited 



molecular electronic states. Since the photoassociation re- 
action occurs at large distances (R > 90 ao), we neglect 
the i?— dependence, using the asymptotic value D^^ de- 
duced from standard long-range calculations j41| . 

In the rotating wave approximation with the instan- 
taneous frequency, the coupled equations for the radial 
wavefunctions 'Fg^ {R, t) in the ground and excited states 
can be written as |2Uj: 



ifi 



d 



/, 



dt 



I, 



T + V{R) + AiR,t) WLfit) 
\WLfit) 'T + V{R)-A{R,t) 

Eq. H17|l corresponds to a rotating-frame transformation, 
with the angle Lij(t), leading to a "frequency-modulated 
frame" |24 | . The coupling between the two channels writes 



(17) 



\'P^{R,t) 



as: 



WLfit)<WLj—^W„ 
TC 



(18) 



Wl = - Jfo^^,^ = -1\^D;- (19) 
2 2 V ceo 



The diagonal terms of the hamiltonian matrix contain the 
kinetic energy operator T, the mean potential V{R): 

Vexc{R) + Vground{R} 



V{R) 



(20) 



and the R- and t-dependent energy difference A{R,t) be- 
tween the potentials dressed by the instantaneous laser 
energy hu){t) (defined in Eg. I16|l : 

2A{R,t) = 2AL{R)-h^ = 2AL{R)-hx{t-tp), (21) 

where 2Al{R) is the R-dependent energy difference be- 
tween the two electronic potentials dressed by the mean 
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laser energy Tilol and which are crossing in R^: 

round 

(R) - hcuL, (22) 
2ALiRL) = 0, 2Al{R ^ oo) ^ 61' (23) 

The instantaneous crossing point Rc{t) between the two 
electronic dressed potentials coupled by the pulse is de- 
fined by: 

A{Rc{t),t)^{)- Rc{tp) = RL (24) 

4.3 Choice of the pulse parameters 

The choice of the present paper is to fix the carrier fre- 
quency /27r (or , equivalcntly, the detuning (5£* ) , the cou- 
pling Wl (or the intensity 7^), and the chirp rate x i^i the 
time domain. In exchange, we shall study pulses of differ- 
ent durations T]^ (or various spectral widths Slo) leading to 
different swept energy ranges 2?i|x|T(7, maxima Wmax — 
^L^^ of the couphng, and energies E^y-^ax ^ WIt^. 

Fixing Wl, instead of E^ax, leads to situations as var- 
ious as possible for either the photoassociation yield or the 
dynamics of the process. Indeed, the pulses which are stud- 
ied induce quite different excitation regimes, due to their 
characteristics: narrow or broad bandwidth 5lo, short or 
long duration tc , small or large chirp <P" inducing various 
ratios ^ and various energies Emax- The object of this 
work is to analyze the properties of the photoassociated 
wavepacket as a function of the pulse. In this aim we shall 
study the photoassociation yield and the radial distribu- 
tion of population transferred in the 0^ surface at the end 
of the pulse, but also the evolution of the wavepackets after 
the pulse (vibration in the 0^ potential and acceleration 
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to the inner region) . An important goal is to gain informa- 
tion about the degree of coherent control in the photoas- 
sociation process: for example we are interested to max- 
imize the population localized at small distances in the 
Og potential (internuclear distances more favorable to ra- 
diative stabilization of the photoassociated molecules, by 
spontaneous or stimulated emission, towards the ground 
potential surface). 

• Detuning. We present here calculations for photoas- 
sociation with different pulses at the same detuning Sf^=2.656 
cm~^, corresponding to a crossing point Rl = 93.7 oq. 
The detuning is chosen such as — Ey^, the bind- 
ing energy of the wo=98 level in the external well of the 
Og" (6s -|- 6^3/2) potential, which means that at t — tp there 

is a resonance condition between the continuum state ly- 
ing at threshold (E=0) and this vibrational level. Taking 
into account that the initial coUisional state considered 
here corresponds to a very low temperature T=54 /iK, 
and to an energy Eq — ksT much smaller than the vi- 
brational energy spacing in the 0^ (6s + 6^3/2) potential 
( "^^'° j^J"^^ ~ 3300), the initial state Eq is in resonance 
with the vq vibrational level. 

• Intensity 1^ and Coupling Wl - Also, we consider the 
same laser intensity II and then the same coupling Wl for 
all the pulses. But the pulses with different durations tl 
correspond to different pulse energies Epuise ~ Iltl- For 
laser excitation with tt polarization between the electronic 
states •^Z'+(6s,6s) and 0^ (6s + 6^3/2), and neglecting the 
_R-variation of the dipole coupling, the intensity is related 
to the coupling Wl by Eq. lfTO|) . giving WL(a.u.)=9.74 x 
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10~^ I ]^{W/cm'^) 0Tj. In the present calculations the which gives, as for two values of the pulse duration 

laser intensity is II — 120 kW cm^^, giving a coupling tc for a given pair (r^, x). In fact, with the relation (|15|l 

Wl— 0-7396 cm^^. By chirping the pulse, the instanta- one can see that the pulse with the longer duration tq 

neous coupling becomes W{t) — WLf{t) < Wmax- corresponds to the higher rate <P" in the frequency domain. 

• Chirp rate x in the time domain. Focussing. The lin- • Temporal width tl of the pulse before chirping. The 

ear chirp parameter x has been designed in order to achieve, Eqs. H25I27|) show that the existence of real values of ^" 

at a time t — tp + Tmb{vo)/2 (where Tmb{vo)/'2= 125 ps and tc requires: 

is half the vibrational period of the level wo=98), the fo- tl < \ — p , (28) 

V \x\ 

cussing of the excited vibrational wavepacket at the in- ^^^^^^^^ ^n upper limit for the value of which can still 
ternal turning point of the vibrational state vo in the 0^- chosen for a fixed value of x- On the other hand, one 
outer well. The chirp parameter necessary to compensate ^^^^^^ ^^^-^ excitation of the 0" continuum 

the dispersion in the vibrational period of the wavepacket ^^^.-^^ photoassociation process. Such a requirement 
is chosen as x = -27rTreAvo) / [TMvo)? , i.e. adjusted to ^^^^^^^ ^ ^^^^^ ^^^^ i^.^i^^ temporal width 

match the revival period T,™ (see Appendix B) of the res- ^^^^^^^ ^^^^ ^ condition for a spectral width of the 

onant level t;o-98. This leads to x=-4.79x IQ-^ ps-2=- ^^^^^ ^^^^^^^ ^^^^ detuning 6f corresponding to the 

0.28 X 10-" au, and ?ix-0.025 cm"! ps'^. vibrational level resonant at t = tp: 



• Chirp rate in the frequency domain. By combin- 4 In 2 



4 In 2 



= n < 5f =^TL> h—^ (29) 

ing the relation (|15|l with Eq. (|14|l one obtains a second l 

degree equation for <P"- ^'^^ ^^'^ detuning (5£*=2.656 cm^^ and the chirp rate x—- 

tn" ,-4 4.79 X IQ-'^ ps"^ leading to focussing, the conditions (|28|l 

and 129() give the following interval of choice for r^: 

This means that for an initial pulse of temporal width 

5.54 ps < TL < 17 ps, (30) 

Tl , there are two values of the chirp <l>" in the frequency 

domain corresponding to the same chirp x in the time corresponding to predominant excitation of bound vibra- 
domain: tional levels 0^ . 



1 / Ixl , /, x'^T'f \ For < Tl < 5.54 ps, dissociation continuum states 

^ ' '-'g localized at large distances are mainly excited. Let 



which give two different temporal widths r^, as it can be 
seen with the formula H14() . • Pulse duration tc. tl and x 
being chosen, is solution of the second order equation: 
, (41n2)2 „ (41n2)2 , ^ 



us remark that the frontier between excitation of bound 
or continuum states is not strict, it only corresponds to 
the FWHM of the spectral energy distribution, without 



T^2^2 c ^ ^2 ^ ' taking into account the absolute intensity of the pulse. 
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In the interval described by the relation (jSOJ, the excita- 
tion of continuum will become important for an increasing 
intensity 1^. The characteristics of the pulses studied in 
the present paper are reported in the Table 1, which also 
shows those of the pulse studied in Ref. |20|. 

4.4 The energy range corresponding to a large 
population transfer induced by the chirped pulse 

The impulsive approximation |30I20| can be used to qual- 
itatively determine the spatial range in which the popu- 
lation transfer from the ground to the excited surface is 
important. Let's take the example of a two- level system 
with energy splitting hueg = h{uje — iVg) which is excited 
by a pulse characterized by a spectral energy distribution 
(centered on the mean laser pulsation lul, FWHM 
Slo, and maximum proportional to the laser intensity I). 
Longtime after the end of the pulse, and in the perturba- 
tive approximation, the population in the excited level is 
given by /igg|£(weg)p, where ^Xge is the dipole moment of 
the transition. Therefore this population grows up when 
fige or I increase and when the detuning Lo^g—LOL decreases, 
being maximum at the resonance ( uj^g = ljl)- 

Similar considerations apply for population transfer 
from the ground to the excited surface; it is favoured: 
(i) at resonance (i.e. at the instantaneous crossing point 
Rc{t)), (ii) for a large overlap integral | <^ Z'+i?o|Og u > ^ 
between the initial coUisional wavefunction and the reso- 
nantly excited vibrational level v having its outer turning 
point close to Rc{t), and (iii) with the increase of the laser 
intensity. 
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For a gaussian pulse, 98 % of the energy is carried dur- 
ing the time interval [~tc, +tc], and then the effect of the 
pulse can be estimated by analyzing it during this tempo- 
ral window |2L)| . During this time interval, the instanta- 
neous crossing point Rc{t) defined by the relations 
varies in time, describing a resonance window] the energy 
range resonantly swept by the pulse around the central 
frequency uj^jlix is 2?i|x|T(7, being limited by its spectral 
width 20 . Indeed, \x\tc < with a ratio depending on 
the chirp rates: 

Then, vibrational levels lying in the "resonance window" 
will be excited by the pulse at different moments. 

Besides, the large overlap of the initial continuum wave- 
function with the vibrational wavefunctions of levels 
close to the dissociation limit is especially favorable to 
population transfer at very large distances. Indeed, the 
overlap integral has a large maximum for vibrational lev- 
els having outer turning points at large distances R, close 
to the dissociation limit (see Appendix C, Fig. I12|l . 

Furthermore, for a large laser intensity, off-resonant 
excitation can easily become efficient. This is particularly 
true for cold photoassociation at small detunings if the 
coupling and detuning have comparable values, leading to 
noticeable excitation outside the resonance window. 

Then, it is interesting to note that using a chirped 
pulse for such free-bound transitions close to the dissocia- 
tion limit, it is possible to control the excited energy range 
by the pulse characteristics, in order to avoid the popu- 
lation of the continuum or to restrain the spatial range 
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Fig. 3. Time-dependent gaussian envelopes f(t) of the linearly chirped pulses considered in this paper, all having the same chirp 
rate x = —4.79 x W~^ps~^ in the time domain, but different tl values, stretched to tc- At each chosen tl value correspond 
two values tc, as explained in the text. The maximum of the pulse, f{tp) = decreases when the time width tc increases, 

a) Chirped pulses with tc « tl= 1 ps, 6 ps and 12 ps, centered at ip=150 ps. b) Chirped pulse with tc=581 ps {tl= 1 ps), 
centered at tp=1500 ps. c) Chirped pulses with rc=96 ps {tl= 6 ps) and rc=47 ps {tl= 12 ps), both centered at tp=1000 ps. 
The horizontal broken line indicates the value from which the pulse becomes efficient in the population transfer. For the same 
intensity II of the laser, the pulses having the same initial tl correspond to the same energy Epuise= "^^t^\/i^i carried by 
the field. 
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covered by the final 0^ packet. Indeed, in addition to the 
spectral width 5uj (determined by for both transform 
limited or chirped pulses) , one can select the range of the 
window IxI'J'c to be excited resonantly, by choosing the 
chirp parameters x, or, equivalently, the ratio ^ (see 
Eq. (El). 

Lastly, the pulse intensity II will determine the adia- 
batic or non-adiabatic character of the transfer (this will 
be detailed in Sec. I4.5|l . and the importance of the off- 
resonance excitation. Nevertheless, the estimated bound- 
ary between resonant and off-resonant excitation remains 
approximative . 

4.5 Condition for an "adiabaticity window" during the 
pulse duration 

Efficient adiabatic population inversion can be obtained 
by using a chirped laser pulse, which in a two state system 
allows to sweep the instantaneous frequency u){t) from far 
above (respectively far below) to far below (respectively 
far above) resonance. Sufficiently slow sweeping induces 
total adiabatic transfer from one state to the other one 

In a previous work |20| we have analyzed adiabatic 
population inversion within the impulsive limit |3U) . as- 
suming that the relative motion of the two nuclei is frozen 
during the laser interaction, i.e. tc ^ Tyif,{vo), where 
Tvib{vo) — 250 ps is the vibrational period of the 0~ level 
resonantly populated at t ^ tp. By neglecting the kinetic 
energy operator appearing in the two-level Hamiltonian 
of the Eq. H17|) , one can introduce a coordinate-dependent 



Cold atoms pliotoassociation using chirped laser pulses 17 

two-level model able to define the conditions for full adi- 
abatic population transfer. At the instantaneous crossing 
points Rc (t) , the nonadiabatic effects can be explored us- 
ing the adiabaticity condition, which in these points takes 
a simple form determined by the pulse shape (|x| and ■^) 
and its intensity {II or Wl) EDI- 

h%\^8{W{t)f. (32) 

The condition l|32|l will not be verified when W(t) be- 
comes very small. Then, it is useful to estimate the do- 
main [—arc, arc] , with a > 0, for which the adiabaticity 
condition (|32|l is verified. For |i — ip| < arc, the coupling 
W(t) has the lower bound -^Wmax, so a can be deduced 
from the condition: 

'-^f^'\xr-^«Wl (33) 

giving: 

" « \/ ln{872^^W^|— } = ama. (34) 
y In 16 h^\x\ TC 

For a <C oimax, the adiabaticity condition H33|l is very 
well satisfied during all the time interval [—arc, arc]. 
For a ~ OLmax, the transfer can be adiabatic during a 
certain portion of this interval, but strong non-adiabatic 
effects appear at the boundary \t — tp\ « otmaxTc- For 
8-f-7^Wl^ < 1 (small intensity II and/or large chirp 
Ixl^), adiabaticity never occurs, amax not being defined. 
The values amax corresponding to the pulses studied here 
are given in the Table 

For a given pulse, amax offers an approximate evalu- 
ation for the extension of the time interval [—arc, arc], 
with a <C amax, during which the adiabaticity condition 
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Table 1. Parameters of gaussian pulses considered in this work, linearly chirped with the same rate x in the time domain, 
and corresponding to the laser intensity at t — tp for the transform limited pulse II = 120 kW cm~^, which, for a linear 
polarization of the electric field, gives a coupling Wl=0.7396 cm~^ for the a'^E^{6s,6s) Oj (6s, 6^3/2) transition at large 
interatomic distances. The central frequency of the pulse, iJL/2TT, resonantly excites the vo = 98 level of the 0~{6s, 6P3/2) outer 
well, with binding energy Evg — 2.656 cm~^ and vibrational period T„i(,(wo) ~ 250 ps. The significant parameters listed below 
are: temporal width tl of the initial pulse before chirping, spectral width hSu, chirp rate in the frequency domain, temporal 
width Tc of the chirped pulse, the ratio the chirp rate x in the time domain, the energy range ?!.|x|''"c resonantly swept by 

the pulse during the period [~rc,rc], the maximum coupling Wmax ~ Wl^J"^ ~ W{tp), and the parameter Qmai indicating 
a limit for the adiabaticity range [—arc, arc] (a <^ amax)- The probability of photoassociation in the Og surface at the end of 
the each pulse (corresponding to a total population normalized at 1 on the grid), noted P. - (Eq), is also shown {Eq is the energy 
of the initial continuum state, Eo/ks = 54 fiK). The probability per pump pulse that a given pair of atoms at the temperature 
T — Eo/kn to be photoassociated is roug hly equal to P^- (T)= P^- (-So) (qe/I rf)|g i ~ 20P(j- (Eq)^, where kpT is the width 
of the thermal distribution and {dn/dE)\EQ the density of collisional states at the energy Eo [Z being the partition function, 
see Eq.Q) ■ The characteristics of the chirped pulse studied in Ref. |20| are reported in the last line. 





hliu) (cm ^) 


<!>" (ps^) TC (ps) ^ 


X (ps-^) 


h\xVc (cm ^) 


Wm.ax (cm ^) 






TL= 1 ps 


14.72 


-0.00063 1.002 0.999 
-208.8 581 0.042 


-4.79X 10"^ 


0.025 
14.52 


0.739 
0.031 


1.12 


3.961 X 10^2 
3.941 X 10"^ 


TL= 6 ps 


2.453 


-0.81 6.012 0.999 
-207.99 96.28 0.249 


-4.79X 10"^ 


0.15 
2.41 


0.739 
0.184 


1.12 
0.5 


2.168 X 10~2 
6.005 X 10^2 


TL= 12 ps 


1.227 


-13.84 12.41 0.983 
-194.96 46.6 0.507 


-4.79X 10"^ 


0.31 
1.17 


0.727 
0.375 


1.11 
0.87 


3.046 X 10"'' 
5.340 X 10"'' 


Tl= 15 ps 


0.981 


-38.51 16.60 0.950 
-170.00 34.8 0.657 


-4.79X 10"^ 


0.42 
0.87 


0.703 
0.486 


1.10 
0.95 


1.200 X 10"" 
3.245 X 10"" 



can be verified at the instantaneous crossing points. Then, 
for fixed |x| and VF^, it is the ratio ^ which fixes the 
time-interval during which crossing of the two potentials 
can lead to population inversion due to adiabatic passage 



5 Results for pulses having different durations 

We consider pulses generated with three initial temporal 
widths: 1 ps, 6 ps and 12 ps, and which arc chirped 
with the same chirp rate in the time domain: x=-4.79x 
10"'^ ps~^. The characteristics of the pulses are given in 
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the Table 1. The impulsive limit is valid for all these 
pulses, except the pulse with rc=581 ps {tl— 1 ps). 

As discussed before, to each given value (determin- 
ing the spectral width 5uj ^ 1 /tl and the total energy of 
the pulse Epuise ^ Iltl) correspond two chirped pulses 
with very different characteristics: one has tc ~ tl, being 
"almost no-chirped" , and the other is a "really chirped" 
one, having a much larger temporal width than the ini- 
tial pulse. The envelopes of the 6 pulses thus built are 
shown in Fig. |21 For each pulse we shall analyze the re- 
sulting photoassociation dynamics: the relative population 
yield transferred to the 0^ state (for a total population 
normalized at 1 on the grid), its radial distribution, and 
the modification of the density probability in the initial 
Eo/kB = 54 fiK a^Z'+(6s, 6s) state. Our aim is to analyze 
their efficiency for the photoassociation reaction, the main 
criterion being a maximum transfer of population from the 
initial (6s,6s) Eq continuum state to bound vibra- 

tional levels of the 0^ (65,6^3/2) excited state: this means 
that we are interested in those pulses able to bring the 
population located at large distances in the initial con- 
tinuum and to put it in the 0^ (65,6^3/2) external well, 
in such a manner that the most of this population to be 
present at small distances (in particular at the inner turn- 
ing point of the outer well) with a small delay after the end 
of the pulse, or to be efficiently accelerated to the inner 
region (in the case of the population transferred at large 
distances). Maximizing the population transferred to the 
inner region of the excited potential (i? < 100 ao) is a first 
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step in a process leading to an efficient formation of cold 
molecules. 

The comparison between the pulses described before 
can be made from two points of view: i) one can compare 
"really- chirped" and "almost no-chirped" pulses, by look- 
ing at results given by the pulses with the same tl (the 
same spectral width hSuj and the same energy Epuise), and 
ii) one can analyze the main differences between the re- 
sults given by short pulses with a large spectral bandwidth 
and those given by much longer pulses with a narrower 
spectral width. 

Several comments are to be made about what we could 
expect from these pulses, given the values illustrated in the 
Table 14.41 and the fact that the detuning corresponding 
to the 0~ (65,6^3/2) level vq = 98 excited ai t — tp is 
5£*=2.656 cm-i. 

First, we can note that the shorter pulses with tl = \ ps 
and Ti=6 ps have the spectral widths hSui— 14.72 cm^-'^, 
much larger than (52*^2.656 cm^^, and 2.453 cm^^, of the 
same order of magnitude as (5£*. This means that, for t^—I 
ps, the 0~(6s,6p3/2) continuum will be massively excited. 
The pulses with 12 ps have the spectral bandwidth 
hSui— 1.227 cm~^ smaller than ^2*, which avoids to notice- 
ably populate the 0^ continuum (except if the intensity is 
too strong, see Sec. I4.4|l . 

Second, the horizontal line in Fig. |31 indicates the value 
from which the instantaneous coupling strength WLf{t) 
becomes efficient in the population transfer. The "efficient 
time durations" of the studied pulses are generally of the 
same order of magnitude as their corresponding tc, ex- 
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cepting the pulse with tc ~ 581 ps {tl — \ ps), for which 
its active interval is reduced by a factor of w 4: indeed, 
due to the strong diminuation of Wmax by the chirp (small 
ratio \J^) only a small fraction of the pulse represented 
in Fig. can produce population transfer. 
Another aspect refers to the value of the detuning 5'^ /2 ~ 
1.328 cm~^ which is always larger than the maximum cou- 
phng Wmax of the studied pulses (see the Table 1). There- 
fore at t = <p , the off-resonant Rabi coupling at very large 
distances i?oo EOI- 

nnitp, R^) = ^wi,,, + {^f (35) 

will be larger than the resonant coupling Wmax at Rl 
(indeed, for the pulses listed in the Table 1, the values 
of fif}{tp, Roo) are between 1.3 and 1.4 cm^^, the corre- 
sponding Rabi periods Tpau = 7r/i7 varying between 11 
and 12.5 ps). This means that, for this small detuning, 
the field is strong enough to couple the two channels at 
very large distances, well beyond we could expect from 
the estimated energy range 2?i|x|rc resonantly swept by 
the pulse during the period [—tc,tc]. In fact, the off- 
resonance excitation at very large distances is enforced by 
several factors: the small detuning S'j^ (corresponding to 
t ~ tp), the coupling Wl, but also the much bigger over- 
lap of the vibrational wavefunctions 0^' v with the initial 
continuum a^S^ (see Appendix C, Fig ll2|) . which for 
u 160 ^ 170 (S,„ = -0.4 -> -0.2 cm-i) is about 
3.7 times larger than the overlap with the state vq = 98 
(i?uf, — —2.6 cm^^). Lastly, the negative chirp presently 
studied begins resonant excitation at t < tp from the 
higher vibrational levels w > 98 of the 0^(6s, 6^3/2) state. 
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These weakly bound vibrational levels, which are excited 
before the maximum of the pulse ai t — tp, are much 
more sensitive at the presence of the field than the lev- 
els excited at t > tp. Indeed, the time evolution of the 
population excited in these states shows Rabi oscillations, 
which are the signature of significant nonadiabatic effects, 
as discussed in Ref. PU]. 

5.1 Population transferred to the 0~ (6s, 6^3/2) excited 
state 

Wc shall analyze the population transferred in the 0^" (6s, 6^3/2 
state during the photoassociation process, using the pulses 
whose envelopes f(t) are represented in Fig.|21 The results 
shown in this section correspond to a total population 
(initially in the a'^Z'^(6s, 6s) state) which is normalized 
at 1 on the whole grid of extension Lp. We represent the 
following quantities: 

• the time evolution of the probability for population 
transfer on the whole grid in the 0~ (6s, 67)3/2) surface: 

P {t)= f'' {R\t)\^dR' (36) 
Jo 

The relative population P^- (Eq) transferred in the 0~ (6s, Gp^/ 
channel after the end of the pulse {t — tp ^ tc) is given 
in the Table 1. The probability of photoassociation per 
pump pulse of a pair of atoms at the temperature T = 
Eo/kp is roughly equal to V,- (T)= P,- {Eo)j^^j§^i 
w 20Pq- (Eq)^, where kpT is the width of the thermal dis- 
tribution and {dn/dE)\Eo the density of coUisional states 
at the energy Eq (Z being the partition function, see 
Eq.Q) 
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• the distribution of the 0^(65,6^3/2) radialy inte- 
grated population, P^- {R,tfoc), as a function of the dis- 
tance R, at the moment t/oc = tp + ^ki. — ip-|-130 ps, 
corresponding to the focussing at the inner turning point. 
At this time moment almost all the studied pulses are 
practically "finished", W{t) ~ 0, except the very large 
pulse with Tc— 581 ps {tl— 1 ps), for which anyway the 
efficient time leading to significant transfer is over (see the 
dashed line in Fig. |3J- For t > tjoc the population P^- re- 
mains constant, the two potentials 0^ and a^^u being no 
longer coupled. 



P,-{R,tf,c)^ \%-{R',tf,,)\HR' (37) 



The Figures 0] and |S1 show the time evolution of the 
population Pt^-{t) transferred in the 0^^(65,6^3/2) state 
during the photoassociation process, using the "almost no- 
chirped" pulses with tc ~ tl = \, 6, 12 ps, and the "really 
chirped" pulses with tq— 581, 96, 47 ps, respectively. On 
the other hand, in Fig. (Hjis represented PQ-{R,tfoc), the 
0~(6s,6p3/2) radialy integrated population as a fonction 
of the distance R, at the moment tjoc ^ tp + -^-^ — 
ip-t-130 ps, for all the pulses. Results for pulses having the 
same tl are represented together, the left column showing 
results until Lji = 19250 ao (the limit of the spatial grid), 
and the right column showing the repartitions of the same 
populations at small distances, until R = 140 ap. 

At a first view the time evolutions of the 0^ population 
given by the three classes of pulses, with r^^l, 6, 12 ps 
(and with three increasing energies of the pulse Epuise) 
are quite different, and this can be understood looking at 
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Fig. 4. Time evolution of the 0g"(6s, 6P3/2) population (for 
a total population normalized at 1 on the grid) obtained by 
photoassociating with the "almost no-chirped" pulses having 
TC ~ TL, whose envelopes f(t), represented in Fig. |3t)i ^.re 
centered at tp — 150 ps. 

the characteristic spectral bandwidths (see the Table 
and comparing them with the detuning. 

As it was mentioned, the pulses with r^— 1 ps will mas- 
sively populate the 0~ continuum (see Figs.QJi) and[Sli)), 
leading at the end to the same large value of the 0^ pop- 
ulation (Pq-=0.04). In this case, the spectral width hSuj 
and the pulse energy Emax appear as the only parame- 
ters playing in the results, independently on the energy 
range 2?i|x|Tc swept during the pulse. The distributions 



22 



E. Luc-Koenig, M. Vatasescu, and F. Masnou-Seeuws: Cold atoms photoassociation using chirped laser pulses 



Cl'^^(R^t +T .72)l^dR' 

^ p vib ' 




R(a^) R(a^ 



Fig. 6. Distribution of the 0~ radialy integrated population P^- {R,tfoc) at tfoc = tp + = tp+130 ps, for "really-chirped" 
and "almost no-chirped" pulses. Results obtained with pulses having the same tl (the same spectral width fiSuj) are represented 
together: thick lines for "really-chirped" pulses (tc S> rz,), and thin lines for "almost no-chirped" pulses (rc ~ tl)- Left column: 
distribution of the Og population until Lj{=19250 ao. Right column: population until R=140 ao. "Really-chirped" pulses give a 
much larger population transfer in the Og state, as the spanned energy range 2?i|x|Tc is much larger; they give also much more 
population at small distances. (The results show the population in the 0" channel relative to a total population normalized at 
1 on the grid.) 
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Fig. 5. Time evolution of the Og population (for a total popu- 
lation normalized at 1 on the grid) during the photoassociation 
process, by using the "really chirped" pulses whose envelopes 
f(t) are represented in Fig. I^tjc): a) Og (6s,6p3/2) population 
for a pulse with tc~ 581 ps {tl~1 ps), centered at tp = 1500 
ps, b) 0^(6s,6p3/2) population for a pulse with tc~ 96 ps 
{tl~G ps), centered at tp = 1000 ps, c) 0^ (63,6^3/2) pop- 
ulation for a pulse with tc= 47 ps {tl~12 ps), centered at 
tp = 1000 ps. 

of the Og population function of the distance R are nearly 
identical for both pulses (see Fig. Et) ), except at short 
distances (Fig.|HJi)). It appears that due to the high value 
of the spectral width hSto > (52', the population is mainly 
transferred at very large distances R > 1000 ao (non reso- 
nantly for the "almost no-chirped" pulse with tc ~ tl = 1 



ps). Only a little amount is excited at smaller distances 
R < 100 ao, where the chirped pulse with tc ~ 581 ps 
has a large "resonance window" (large value fi.|x|Tc) and 
it is more efficient than the short pulse of rc ~ 1 ps, with 
a narrow "resonance window" . It is interesting to remark 
that, for the pulse with rc w 1 ps the population transfer 
is adiabatic {amax > 1) during the time window [~tc, tq], 
but the adiabaticity condition cannot be satisfied for the 
very long pulse with tc — 581 ps, for which amax is not 
defined (see Table 1). In fact, the oscillations which can 
be observed in the time evolution of the 0^ population 
in Fig. [SJl) correspond to a strong Rabi coupling at very 
large distances (described by the Eq. (|35|l ) and are the 
signature of a strong nonadiabatic behaviour in the pop- 
ulation transfer. This non adiabatic population transfer 
at the large distances swept by the instantaneous crossing 
point results in a large population remaining in the 
surface after the end of the pulse. Similar oscillations in- 
dicating a nonadiabatic transfer can be equally observed 
in Fig.jSJs), in the evolution of the 0~ population during 
the excitation with the pulse with tc = 96 ps {tl— 6 ps), 
to whom corresponds a small value of Umax = 0.5. 

The two other classes of pulses, with tl= 6 and 12 ps, 
have much narrow spectral widths, smaller than the de- 
tuning (5£*=2.656 cm^^. For these cases, it appears clearly 
that, for the same r^,, the "really chirped" pulse is much 
more efficient for the Og" population transfer, both for the 
total transfer, as for the transfer at small distances (see 
Figs.|nb,c,e,f). The radially integrated population P^- {R,tfo 
is generally a linearly increasing function of the distance R 
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(see Fig.Elleft column) except for the pulse with rc= 47 
ps {tl— 12 ps), for which Pg- {R,tfoc) reaches its limiting 
value at small R w 140 ao (Fig-Et,f)), because the popu- 
lation transfer takes place adiabatically in a "photoassoci- 
ation window" [201 • Indeed, in this case, the energy range 
resonantly swept by the pulse to the large distances during 
the period [—tc^tc] is TiIxItc ~ 1-17 cm^^, the spectral 
width is hSuj =1.23 cm^-*^ (satisfying ?i|x|rc + fiSu; < (5£*) 
and the population transfer keeps an adiabatic character 
{oimax — 0.87) both in the resonance window, and at large 
internuclear distances (^(-R) ^ Sfj^ for R > 200 ao). 

On the contrary, the chirped pulse with tc= 96 ps 
{tl= 6 ps) does not produce excitation in a photoasso- 
ciation window, but everywhere at large distances, be- 
cause the energy range swept during the time window 
(to the large distances) approaches the dissociation limit 
{h\x\Tc = 2.41 cni-i w 61*). 

We have to remark that the very different evolutions 
of the 0^ populations during the pulses with tl= 6 ps 
(see Figs. 0Jd), |S|d)), on the one hand, and with tl= 12 
ps (see Figs. 01;), Et)) on the other hand, are typically for 
two kinds of results that can be obtained: at the end of the 
pulse one can have either population excited everywhere 
at large distances (in an amount increasing with R, re- 
flecting the amplitude of the initial contimmm), or, if the 
spectral band of the pulse is narrow enough, one can ob- 
serve a selective adiabatic transfer of population in a cer- 
tain domain of distances, representing a "photoassociation 
window" , and very small population transfer outside the 
window (where the transfer is also adiabatic, but where 
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Fig. 7. Excited wavepackets localized at tfoc = tp + 
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at the inner turning point of the outer well of the 0^ poten- 
tial. Comparison of two typical final states: a) transfer at large 
distances, and b) excitation in a "photoassociation window". 



there are no instantaneous crossings, which explains the 
small value of the remaining P^- at the end of the pulse). 
An illustration of these results is made in Fig. [3 showing 
the Og wavepackets at tfoc = tp + = tp+ 130 ps, a) 
for the chirped pulse with tc~ 96 ps {tl~ 6 ps), giving 
population transfer on the whole spatial grid, and b) for 
the chirped pulse with tc— 47 ps {ti^— 12 ps), producing 
transfer in a "photoassociation window" . Due to the choice 
of the X value, at tjoc = tp + both wavepackets 

are focussed at the inner turning point of the vibrational 
eigenstate vq excited at t = tp. 
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We are especially interested in pulses producing the 
maximum transfer of 0^ population to the small and in- 
termediate distances (see Fig|Sl for this analysis). From 
this point of view, the pulses with tl= 6 ps {tc— 96 
ps and Tc— 6.012 ps) and the pulse with tc= 47 ps 
{tl— 12 ps) are giving the bigger population at small R. 
Their results can be seen and compared in Figs. Et) if)- 
tfoc = tp + both pulses with tc— 6.012 ps {tl= 6 
ps) and with tc— 47 ps {tl— 12 ps) give Pq- (100 ao, tfoc) 
« 4 X 10^"*, but the pulse with tq— 47 ps, which produces 
a "photoassociation window" , is more favorable to transfer 
population at smaller R-values R < 100 ao- Nevertheless, 
it is the pulse with tc— 96 ps {t^— 6 ps) which seems 
to be the most efficient for the 0~ population transfer, at 
large distances, but also to the inner region (Pg- (100 ao, 
tfoc) ~ 8 X 10~^), because even after the end of the pulse, 
the population will be accelerated inside from the large 
distances. 

Let's remark that, as dicussed in Ref.|2ni) the negative 
chirp has a noticeable contribution in the acceleration of 
the wavepacket created in the excited state 0~ towards 
short intcrnuclcar distances. Indeed, Lu{t) decreasing with 
time, the instantaneous crossing point moves to smaller 
distances, following the motion of the wavepacket which 
is accelerated inside the potential well. This is favorable to 
our goal of maximizing 0^ population at small distances. 
On the other hand, this means that an interpretation of 
the transfer process within the impulsive approximation 
could easily be invalidated. For example, in our case, the 
population in the initial continuum state integrated until 
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R = 100 ao is: 

I-R=1Q0 

/ \^s.T=54^.K{R')\''dR' = 3 X 10-4 (38) 
Jo 

and this is what would be the maximum value of the in- 
tegral /o^~^°" 1^0^ (-^'' tfoc)\^dR' if the impulsive approxi- 
mation would be valid until tfoc- But we have shown that, 
in the cases of some of the pulses discussed before, (with 
Tc «6 ps, Tc= 96 ps and tc= 47 ps), Pq-(100 ao, tfoc) 
reaches indeed bigger values, proving the acceleration to- 
wards shorter internuclear distances for times t < tfoc- 

5.2 Population in the last vibrational levels of the 
ground state 

In the non perturbative regime there is a redistribution of 
population in the ground surface a^Z'+(6s, 6s), arising as 
well in the bound spectrum, as in the dissociation contin- 
uum. There appears a "hole" in the initial wavefunction 
around Pl, a large part of this population being trans- 
ferred in bound vibrational levels of 0^ and a^S^ - 

Due to the coupling at large distances between the 
two electronic states, the last vibrational levels of the 
a^IJ^(6s,6s) state are noticeable populated during the 
photoassociation process. Figs. |H1 and |51 show the evolu- 
tion of the population in these last levels {v" = 53 is the 
last vibrational level [20]). If one compares pulses with 
the same r^, the final transferred population is roughly 
of the same order, but for tl= 1 or 12 ps it is bigger for 
"almost no-chirped" pulses, probably because the maxi- 
mum coupling Wmax is more intense for shorter pulses. 
The pulses with tl— 6 ps {tc~ 6.012 ps and 96 ps) give 
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almost the same population of v"=52, 53. In Fig.lHh)) we 
show the Rabi oscillations in phase opposition between the 
Og population P^- (which is distributed on the whole spa- 
tial grid of ii^=19250 ao) and the population of the last 
vibrational state of the a'^Z'^(6s, 6s), with v"= 53, whose 
wavefunction extends at very large distances (the last os- 
cillation is between 150 and 1200 ao, see Fig. |2K))- This 
result is symptomatic for the efficient exchange of popu- 
lation between the two channels at very large distances, 
already emphasized, and due to the small detuning and 
large overlap integral (see Appendix C, Fig. I12|l . Around 
tp — 1000 ps, the period of the observed Rabi oscillations 
is about 12.4 ps, in agreement with the result given by 
the formula l|35(l . The oscillations which can be observed 
in Figs. and b) are characteristic for nonadiabatic ef- 
fects in the transfer of populations. 

One has to emphasize that, in the cases of adiabatic 
transfer in a limited spatial window (exemplified here with 
the pulse tq— 47 ps, and in Ref. 20' with a pulse hav- 
ing Tc~ 34.8 ps), the populations in the last bound states 
v"= 52,53 of the ground potential at the end of the pulse 
is of the same order as the population P„- remained in 
the G~ levels resonantly excited by the pulse in the so- 
called "photoassociation window". For example, for the 
pulse with Tc= 47 ps, P^- (Eq) ^ 5.34 x IQ-* (gee Ta- 
ble 1), and the population P^s^iv" = 52) -I- ft^^_(u" = 
53) w 2.25 X 10-4 (see Fig.^t) ). Then, in this case, im- 
mediately after the pulse, the individual 0^ bound lev- 
els belonging to the "photoassociation window" are less 
populated than these last a'^Z'+(6s, 6s) states. But it has 
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Fig. 8. Populations in the last vibrational levels 
v" (a'^X'^ (6s, 6s)) of the ground state, for "almost no-chirped" 
pulses (rc ~ tl). (The total population is normalized at 1 on 
the grid.) 



to be noted that v" =52,53 have wavefunctions extend- 
ing at very large distances (hundreds of ao), being then 
well populated during the exchange of population between 
the two channels, favoured by the good overlap at large 
distances, as explained before. In exchange, the levels 0" 
populated in a "photoassociation window" are bound at 
much smaller distances, as their outer turning points are 
between i?,„i„=84 ao and ao. 
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Fig. 9. Populations in the last vibrational levels 
v"(a'^£'^ (6s, 6s)) of the ground state for chirped pulses 
with Tc= 581 ps, 96 ps, and 47 ps. (The total population is 
normalized at 1 on the grid.) 

5.3 Evolution of the 0^ wavepacket after the pulse 

We shall briefly describe the wavepackets evolution in the 
Qg surface after the pulse, in the two typical cases dis- 
cussed before: excitation in a limited spatial "photoasso- 
ciation window" and excitation at all distances, with mas- 
sive transfer of population at large distances. Obviously, 
these two types of wavepackets created by the photoasso- 
ciating pulse are extremely different, and their dynamics is 
significant in view of some anticipation of the results that 
could be brought by a second pulse. This second pulse can 
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Fig. 10. Og wavepackets evolution after the pulse with rc = 47 
ps {tl = 12 ps). 

be even identical with the first one, but time-delayed (in 
the analysis of the repetition rate of the laser) or a differ- 
ent pulse, if the goal is the stabilization of the system by 
stimulated emission to low vibrational levels of the ground 
state. These subjects will be treated in a future article. 

5.3.1 Excitation in "a window": Vibrational dynamics. 

Fig-Elshows the evolution of the 0^ wavepacket adiabat- 
ically excited {oimax — 0.87) in a limited spatial window, 
with a chirped pulse of time width tc — 47 ps (r^ = 12 
ps). In fact, during the pulse, a large amount of population 
is transferred to levels close to the Cs2(6s-f'6p3/2) disso- 
ciation limit (see the wavepacket at t = <p, in Fig. llOh ) 
), but, due to the adiabaticity of the population transfer 
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outside the "photoassociation window" , this population 
goes back to the ground state and only a small range of 
distances (R<200 ao) remains populated. Then, after the 
pulse, one can observe a typical vibrational dynamics, as 
it can be seen in the figure for some selected moments: 
at t = tp + '^"'''^'■""^ the packet is focalized at the inner 
turning point of vo, a,t t — tp + Tyib{vo) the packet comes 
back at the external turning point of vq {Rl = 93.7 oq), at 
t = tp + 2Tyii,(vo) ~ tp + 500 ps, one can distinguish two 
parts: one, with a maximum around = 93.7 oq, is com- 
posed by levels with vibrational periods close to Tmi,{vo)', 
the other has a maximum around R = 130 ao- indeed, the 
vibrational levels of the 0^ potential having the external 
turning point around this distance vibrate with Tyih ~ 500 
ps. 

5.3.2 Excitation at large distances: acceleration to the inner 
region 

For pulses having sufRciently large bandwidths {Suj > (^2'), 
the population transferred at large distances during the 
pulse remains on the excited surface 0~ also after the end 
of the pulse. Just after the end of the pulse, the radial 
distribution j'Z'g- (i?, <)p of the probability density repro- 
duces that of the initial coUisional state on the sur- 
face at sufficiently large distances R > 500 ap. After that, 
the wavepacket evolves in the —C^/R^ potential and it 
is accelerated toward the inner region. Fig. 1111 describes 
the evolution of the wavepacket created by the pulse with 
Tc = 96 ps {tl — 6 ps), at different moments tp + At, 
with At =0.1, 4, 9, and 14 ns. As it can be seen in Fig. 
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lllb). for R > 1500 ao, there is no noticeable modification 
of the I'Z'g- {R, t)\, even for the largest value of At. As At 
increases, the nodes occuring at i? > 500 ao (i?i ~ 800 
ao and i?2 ~ 1270 ao, for example) begin to be shifted 
toward small R values. In the range R < 500 ao, vibra- 
tional motion can be observed winning progressively dis- 
tances Rvib{tv) in the outer well of the 0^ potential, after 
a time t„ approaching the vibrational period T^n, w t^ 
of the levels whose wavefunctions have the outer turn- 
ing points around Rvib{Tvib)- For example, for v=162, 
with Tyib w 13.7 ns, Rmb = 490 ao. (The largest time 
Tmax — 2.9 /is relevant for the presently used grid with 
Lp = 19250 ao correspond to vibrational movement at 
distances R w 4000 ao.) 

The acceleration of the 0^" wavepacket can be esti- 
mated introducing the radial flux |46|: 

J^-{R,t) = -'hm[^* {R^t)^^^-{RM (39) 

" jJL 9 OH 

with /i = 121136 a. u. the reduced mass of Cs2. Writ- 

i<p _ (R.t) 

ing the wavepacket as ^Q-{R,t) = (i?, t)|e , 
one has \Po-{R.t) = = ^^^_{^R^t). Fig. 

Illb ) reports the R-variation of the momentum ^Po^ 
the times tp + At; the rapid increase at intermediate R 
values 200 < R < 700 ao, for 0.1 ns < At < ^ ns is 
clearly illustrated. This increase in the momentum trans- 
fer during the pulse is a signature of the "kick" given by 
the laser light to the molecule. For "really chirped" pulses 
the photoassociation probability is relatively important, 
and the population transferred at intermediate distance 
is accelerated to the inner region a short time after the 
end of the pulse, increasing significantly the population at 
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Fig. 11. a. 1-4) Og wavepackets evolution after the pulse with tc = 96 ps {tl = 6 ps), at different moments after the end of the 
pulse (whose maximum is at tp = 1000 ps): t= tp+0.1 ns, tp+ 4 ns, tp+9 ns, tp + 14 ns. b) The R-variation of the momentum 
j-Pg-{R,t) at the same times: tp+0.1 ns (dot-dashed line), tp+ 4 ns (thick continuum line), tp+9 ns (dashed line), tp+14 
ns (thin continuum line), showing the acceleration of the wavepacket to the inner region. The singularities in the momentum 
behaviour are due to the definition modulo tt of the phase. 
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short distances. Optimization of the pulse parameters in 6.1 Analytical thermal average using the scaling law in 

order to give to the population transferred at intermediate the vicinity of E ^ 

distances a strong momentum to the inner region will be 



considered in the future. 



In Ref. |2(J| and in the most part of this paper we present 
results of photoassociation calculations considering only 
the s-wave and taking as the initial state in the a^U^{6s, 6s) 
g Yhe Photoassociation probability from a ground potential a single continuum state of the thermal 



thermal average over the incident kinetic 
energies 

Supposing thermal equilibrium at the temperature T, the 
initial density matrix can be expressed in terms of energy- 
normalized coUisional eigenstates of energies E in the ground 
state a^^uj the formula Then the total prob- 

ability per pump pulse that a given pair of atoms to be 
photoassociated into the excited state 0^ at the tempera- 
ture T, is: 



energy distribution, with the energy Eq. In our calculation 
Eq = 1.7211 X 10^1" a.u.= 3.778 x 10"^ cm^^, being, in 
our method, the energy of a continuum state belonging to 
the discretized continuum which is calculated numerically 
with the Sine Mapped Grid method in a box of radius 
Lii = 19250 slq. Then, with this choice of the initial con- 
tinuum state, the temperature in the thermal ensemble of 
atoms is arbitrarily defined hy ksT = Eq, as being T = 54 
fiK. At this temperature, the initial Boltzmann distribu- 
tion (see the formula © ) is narrow (e~^ = 0.05), and a 
first approximation is to represent it by its mean energy 



jy ^j.-^ _ _}_ / dEe~^^ P (E) (40) ^'^^ which is justified by the fact that the spectral widths 

Os Z Jq °9 

hSuj of the pulses considered in photoassociation (see Ta- 
where P,- {E) accounts for the density of probabihty (prob- ^^^^ j^^.^^^. ^^^^ ^j^^^^ ^j^^ ^^^^^^^^ 

ability per unit energy range) that an energy normalized evaluated from P„- (i?o), by using a threshold scal- 

continuum state of energy E belonging to the ground state ^^^.^^^^ probabilities P^- {E) for E ^ E^ 



a S+ to be photoassociated at the end of the pulse in the 



excited state 0„ . Pn-(E) has dimensions of 1/E, being i n j- j-i, j-i- i • tj 

ff 0, V / / 5 6 shall discuss the conditions making valid the use 

obtained with the formula 15()ll from the dimcnsionlcss r;i r i r ^- ^- f^u i- 

' — ' of the scaling law tor estimations of the photc 



toassocia- 



photoassociation probability P„- (i?) in the 0„ state: ^ tit u • j-i j- j- • i 

^ ^ Og y ■> g tion rates. We emphasize that such a discussion makes 

Bjt; sense only for the photoassociation with a cw-laser or 

""o.i^^-io^r'^i^)^ (41) 

with a pulse leading to the formation of a spatially lo- 

We shall discuss in the following the evaluation of the calized wavepacket on the excited surface |l(3l20) . giving 

integral 140|l . what we have called a "photoassociation window" . In such 
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a case, the main contribution to the photoassociation pro- 
cess towards high excited vibrational levels is provided 
by the very localized range of the internuclear distances 
[Rmin , Rmax] i around R ^ Rl , swept by the instantaneous 
crossing point Rc{t) during the time window [— Tc, Tc] f2I)j. 
The photoassociation yield is determined by the overlap 
between the initial stationary continuum P-S"^, -E > in the 
ground state and the excited wavefunctions |0g ,u > (see 
Appendix C, Fig. IH]! . 

A scaling law for the behaviour of the continuum wave- 
functions can be obtained from their asymptotic forms. 
For a very small detuning 5^*, the crossing point i?L of 
the dressed potentials is at a distance large enough mak- 
ing that, even at a low continuum energy E, the potential 
in Rl can be considered as negligibile: E > Then, for 
a s-wave, the continuum wavefunction can be described 
by its asymptotic behaviour |43) : 

\%AR)\ - {^^^^^^^ (42) 

L is the scattering length of the ground surface. This 
means that, for sufficiently low collision energies: 

= (43) 

one obtains the following probability density in Rl'. 

Wg.E{RL)?-%HRL-Lf-^VE (44) 

as E — (hk)'^ /2^. Then, for pulses producing a spatially 
localized wavepacket on the excited surface, and for suf- 
ficiently small energies, one can assume the following law 
for the probabilities pQ-iE) ^ |!f'g,£;(i?L)p: 

P„-{E)=P,-{Eo)^/E/E~o (45) 
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Using the relation (|45|) , the integral (|40|l can be evaluated 
as being: 

P^!}{T) = lp^^iEo)-^ikBT)'/' (46) 

ZV°'-{T) represents the probability that a pair of atoms, 
belonging to a gas in a volume V and at the temperature 
T, and described only by s-waves, to be photoassociated 
in the 0^ (6s -I- 6^3/2) state. 

For a typical number which can be obtained for this 
quantity, we shall take as example the pulse studied in 
great detail in Ref-j^Ol) whose characteristics are given 
in the Table 1, and having the duration tc — 34.8 ps 
(tl = 15 ps). The induced dynamics (which is very simi- 
lar to that resulting from the "really chirped" pulse with 
Tc — 46.6 ps and tl — 12 ps) results in an adiabatic 
population transfer in a spatial range with Rmin = 85 ao, 
Rmax = 110 ao, with amax = 0.95, and populating about 
15 vibrational levels in the vicinity of vq = 98 level in 
the external well of the 0g"(6s + 6^3/2) potential. It cor- 
responds to the photoassociation probability P^- (Eq) — 
3.245 X 10^"*. Then, taking into account the density of 
states j^\E„ = 1.1415 x 10" a.u. at Eo/ke = 54.35 ^K, 
for a gas of cesium atoms at T = 54 /iK we have used the 
formula (|46|l to obtain: 

ZV^-(T = 54 fiK, Tc = 34.8 ps) = 0.00560 (47) 

with Tc characterizing the pulse used for photoassociation. 
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6.2 Average implicitely accounting for real threshold 
effects 

As discussed in a previous section, in the conditions of 
temperature and detuning discussed in the present paper, 
the asymptotic behaviour of the continuum wavefunction 
having the energy Eo — ksT, T — 54 fiK, is not reached 
at i? = Rl = 94 ao. Indeed, the potential energy strongly 
determines the structure of the initial wavefunction for 
R 1^ Rn — 82.3 ao, Rn being the position of the last com- 
mon node. This observation is in agreement with the dis- 
cussion of Ref . j4H| . showing that the asymtotic behaviour 
in a potential is reached for i? > i?^ = (f^)^^^, 
giving Rb — 95 ag for the ■^Z'+(6s,6s) potential. 

Then, in the present studied example of cold atoms 
photoassociation, the detuning 61* = 2.656 cm^^ is too 
large {Rl ~ Rn), and the temperature T = 54 jiK not 
sufficiently small compared to Tan — 69.2 jiK (see the 
relation (|43|l ), to allow the evaluation of the photoas- 
sociation rate from an analytical thermal average. The 
scaling law in \/E being not valid, it is necessary to ex- 
plicitely study the energy variation of P^- (E), by consid- 
ering different initial collisional states of energy E in the 
^Z"^ potential, and having a node at Ljj. In this case, 
the real threshold effects are completely and correctly ac- 
counted for implicitely by performing numerical integra- 
tion in Eq. (|40|l . For example, for the pulse with tc — 34.8 
ps {tl — 15 ps) studied in Ref.|20|. we have performed 
21 time-propagation calculations for collisional energies in 
the range 36.6 nK < E/ks < 633 fiK. Let us emphasize 
that this detailed analysis of the threshold effects has be- 
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come possible owing to the Mapped Sine Grid method j^Sj 
for which a large box of dimension Ln can be considered. 
The energy-dependence of Pq- (E) will be analyzed in a 
further publication . It differs strongly from the 
scalling law. Pp- {E) increases very rapidly at threshold 
exhibiting a very sharp asymetrical resonance like struc- 
ture with a maximum at E/ks ~ 7.9 fiK (smaller than the 
temperature studied presently) with a FWHM AE/ks ~ 
4.7 /iif . For the temperature T — 5A /iK, numerical inte- 
gration in Eq. (|in|l leads to 

ZV^- [T = 54 ^lK, TC = 34.8 ps) = 0.00685 (48) 

The rather good agreement between the analytic average 
(|47|l calculated with the scaling law, and the average (|48|l 
containing the threshold effects is completely fortuitous. 

The two pulses described by the same parameters (52*, 
Wl and and differing only by their durations tc = 34.8 
ps {tl — 15 ps) and tc = 46.6 ps [tl = 12 ps) are as- 
sociated with rather similar "photoassociation windows" 
characterized by the values [amax = 0.95, Tilxkc = 0-87 
cm^^) and {a,nax — 0.87, ?i|x|rc ~ 1.17 cm^^). For the 
same initial collisional state Eo/ks = 54.35 /iK, the cor- 
responding photoassociation probabilities are Pg- (Eq) ~ 
3.245 X 10"* and 5.340 X 10 The photoassociation prob- 
ability is larger for the longer pulse with tc ~ 46.6 ps 
(tl = 12 ps), which creates a slightly larger photoassoci- 
ation window {Rmin — 84 ao, Rmax — 117 ao), compared 
to {Rmin = 85 ao, Rmax = HO ao) for the pulse with 
Tc — 34.8 ps {tl — 15 ps). This increase of Pq- (-Eq) has to 
be related to the increase of the probability density within 
the photoassociation window in the initial collisional state 



E. Luc-Koenig, M. Vatasescu, and F. Masnou-Seeuws; 

(/i!rr I'^s^.EoiRll^dR' = 0.00044 and 0.00076, see also 
Ref. j2()| . Fig. 3), which is a signature of the nearly total 
population transfer in this range of R. 

For the pulse with tc = 34.8 ps (t^ = 15 ps) we 
have shown that the energy dependence of Pq- (E) is pro- 
portional to that of the square of the overlap integral 
\ < Og vn ^ 98\^S+ E > \^ Taking into account 
the very similar dynamics induced by the two pulses with 
Tc — 34.8 ps {tl — 15 ps) and tc — 46.6 ps {tl = 12 
ps), it is reasonable to assume that the energy variation 
of Pq- (E) is the same for both pulses. Therefore, for the 
same temperature, both pulses correspond to the same 
value of the ratio T'p- (T)/Pq- (£'o), and the total prob- 
ability per pump pulse tc — 46.6 ps (r^ — 12 ps) that 
a given pair of atoms to be photoassociated into the 0^" 
state at a temperature T = 54 fiK can be estimated to: 

ZVf,- (T = 54 fj.K, TC = 46.6 ps) = 0.0113. (49) 

6.3 Total number of molecules photoassociated per 
pump pulse 

For a number of N atoms in a volume V, the number of 
pairs of atoms is ^-^-^ — — ~ Taking into account the 
spin degeneracy of the Cs{6^S) atomic state, dA = 2, and 
of the initial electronic state dsi;,, — 3, the total number 
of molecules photoassociated in the excited state 0~ per 
pump pulse is: 

^f--Vo,iT)^ (50) 

For a trap of volume V — 10~^ cm'^, at the temperature 
T — 54: jiK and with a density of atoms Na = N/V = 
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10" cm-3, the partition function is Q{T) ^ 59.86 x 10"^° 
ag"'^^ 40.4 X 10^^ cm^'^. Then the number of molecules 
photoassociated per pump pulse is: 

A^(rc=34.8ps) = 0.69, AA(rc=46.6ps) = 1-40 (51) 

For a repetition rate equal to 10^ Hz and supposing that 
each pulse acts on the same initial state, this gives 6.9x10^ 
molecules per second for the pulse with tc ~ 34.8 ps (t^ = 
15 ps) and 1.4 x 10® molecules per second for the pulse 
with Tc — 46.6 ps (r^ = 12 ps). 

The analysis of the energy-dependence of the photoas- 
sociation probability P^- (E) for pulses leading to signifi- 
cant population transfer at large internuclear distances is 
in progress. 

7 Discussion : possible ways for optimization 

From the previous analysis, we may extract some direc- 
tions on possible ways of optimizing the pulse. In the situ- 
ation where the population transfer occurs mainly within 
the photoassociation window, for sufficiently large cou- 
pling Wmax the total adiabatic population transfer im- 
plies that the whole population for pair of atoms with 
relative distance lying in the [Rmin-,Rmax] range, Pinu = 
jfl^ax |^^^^^(^/^ I ^ 0)|2|di?', is transferred to bound lev- 
els of the excited state and to the last bound levels of 
the ground state. The photoassociation yield can be opti- 
mized by designing the chirped pulse in order to maximize 
Pinit- This can be achieved by increasing the photoasso- 
ciation window with tc values as large as possible, under 
the condition filSu) + |x|tc) < ^2* so that only bound 
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levels are populated. This yields an optimal value for 
which is Topt ^ ^^i? ^ . Since the condition (28) is fixing an 
upper value for to make focalization possible, we end 
with an upper limit for the detuning, and optimization 
can be achieved making use of the scaling laws govern- 
ing the spectra of long range molecules. This will be fur- 
ther explored in future work, but for the detuning (52*= 
2.652 cm~^ considered in the present work, it is clear that 
the optimal pulse would be tl ^ 10 ps, corresponding 
to Pinit^ 1.27 X 10^"^, i.e. to a photoassociation rate in- 
creased by a factor 2.4 compared to the pulse with tc= 
47 ps {tl—12 ps), considered in the present work. 
We have also obtained an increase of the population by 
a factor of 4.3 when considering the pulse rc=34.8 ps 
[TL — \b ps), and increasing the coupling by a factor of 
16, corresponding to a peak intensity 7^=3.36 MW cm"^ 
for the transform limited pulse. 

Finally, a factor of 5 on the probability reported in the 
last column of Table I was obtained in case of the pulse 
with Tc— 47 ps (tl = 12 ps), by reducing the detuning to 
0.0695 cm~^, therefore moving the value of Rl to 150 ao. 

8 Conclusion 

We have investigated the possibilities offered by chirped 
laser pulses to optimize the yield of the photoassociation 
process. Following a previous paper |20j . time-dependent 
calculations have been presented for the particular exam- 
ple of the reaction Cs(6s) -t- Cs(6s) Cs2 0^ (6s 
+ 6P3/2) involving ground state cesium atoms at a tem- 
perature T ^ dAjiK, colliding in presence of laser pulses 
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of different spectral widths Slu with the same linear chirp 
rate in the time domain x=-4.79 x 10 ps^^ . The cen- 
tral frequency is red-detuned by Sl ~ 2.65 cm~^ relative 
to the D2 atomic line, so that the central frequency of the 
pulse excites at resonance the wo=98 vibrational level in 
the outer well of the 0~ level potential, mainly through 
a vertical transition at the distance Rl ~ 93.7 oq cor- 
responding to the outer classical turning point of the vq 
level. The new aspects in the present work are : 

— i) the calculations take into account the mixed state 
character of the initial collision state, described by a 
statistical mixture of stationary collision eigenstates 
representing thermal equilibrium at T '-^ bAfiK . This 
choice reproduces correctly the spatial delocalization 
of the initial state, with large de Broglie wavelength 
^DB ^ 975 ao. In the range of detunings considered 
here, experiments with a cw laser have demonstrated 
minima in the photoassociation rate corresponding to 
the nodes in the scattering wavefunctions : such nodes, 
common to the various eigenstates, are correctly repro- 
duced in the present treatment. The continuum wave- 
functions are represented as eigenstates in a large box 
(size Lj^ ~19 500 ao), using a mapped sine grid repre- 
sentation which involves a reasonable number of grid 
points. The time-dependent Schrodinger equation de- 
scribing motion in two realistic potential curves cou- 
pled by the laser field is solved numerically. 
At ultracold temperatures, the threshold effects which 
govern the dependence of the photoassociation rate as 
a function of the energy E have to be accounted for 
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correctly. A high resolution analysis of such effects can 
be achieved by considering a large number of unity- 
normalized wavefunctions in the box, from which en- 
ergy normalized functions are deduced. A proper esti- 
mate of the absolute value of the photoassociation rate 
has then been obtained from an incoherent average 
over a thermal distribution of the energy normalized 
wavefunctions . 
— ii) In the present paper, we have explored a large va- 
riety of pulses, all of them having the same central 
frequency, as described above, resonant with the level 
vo=98. All of them have the same linear chirp param- 
eter in the time domain, x, designed so that after the 
pulse, at time tp + Tyih/2 (where tp corresponds to 
the maximum of the pulse, while Tyib is the classical 
vibrational period of vq) the vibrational wavepacket 
created in the excited state is focussing at the inner 
turning point. This choice is dictated by the objective 
of improving the efhciency of the stabilization step, 
where either by spontaneous or by induced emission 
the population is transferred to bound levels of the 
a^i7+ ground triplet state. This chirp rate can easily be 
deduced from the revival period, which can be defined 
provided the levels populated stay in a small energy 
range around vq. All the pulses are obtained from a 
gaussian transform-limited pulse, with the same peak 
intensity II— 120 kW/cm^. They differ by the spec- 
tral width Suj ~ (tl)^^. By varying tl from 1 to 6 and 
12 ps, various situations are analyzed, with narrow (~ 
1 cm^^) or broad (~ 15 cm~^) spectral width, long 



Cold atoms photoassociation using chirped laser pulses 35 

(t(7 ^ 580 ps) or short {tc ~ 1 ps ) duration of the 
pulse stretched by chirping. Indeed, for a given choice 
of Slo and x, two different pulses can be associated : 
the first one, with {tc/tl) « 1 is referred to as "al- 
most no-chirped", the second one, with [tc/tl) ^ 1 
as "really chirped" . The resonance window 2?i|x|Tc ex- 
plored by the instantaneous frequency during the time 
interval [tp — Tc.tp + tq] is larger in the second case. 
Different situations are then encountered with respect 
to the validity of the impulsive approximation, or of 
the adiabaticity of the population inversion. 

The conclusion of our study is that the "really chirped" 
pulses seem generally to be more efficient. In fact, we have 
observed two qualitatively different dynamical situations, 
depending upon the energy range defined by h,{5uj-\'\x\Tc), 
i.e of the spectral width and the width of the resonance 
window : 

— Under the condition h{5u) -f Ixkc) < (^J', already con- 
sidered in previous work 20 , the concept of a photoas- 
sociation window is relevant. In spite of the delocalized 
character of the initial wavefunction, we observe after 
the pulse a wavepacket in the excited state localized 
around Rl with a finite extension [Rmim Rmax], corre- 
sponding to the domain of variation of the outer turn- 
ing points of the vibrational levels in the photoassocia- 
tion window. Indeed, whereas during the time window 
[tp — Tc, tp + Tc] levels outside this window are signifi- 
cantly populated, it is possible to optimize the param- 
eter of the pulse so that an adiabatic model is valid 
where no population remains, after the pulse, outside 
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the photoassociation window. By varying the size of 
the photoassociation window, it is possible to optimize 
the number of photoassociated molecules, and this di- 
rection should be further explored in future work. 
For fixed detuning and spectral width, an increase of 
the laser coupling Wl modifies the photoassociation 
dynamics. Having defined a parameter a such that adi- 
abatic population transfer is taking place in the win- 
dow [—arc, +aTc], we have shown that a larger inten- 
sity is increasing the time window, and therefore the 
width of the photoassociation window [— ?iQ!|x|rc, -|-?iQ!|x| 
where total population inversion is taking place. There- 
fore as expected the photoassociation rate should be 
increased at large intensities. However, the dynamics 
outside the photoassociation window at large internu- 
clear distances {R > 200 ao) becomes less adiabatic. 
Rabi oscillations appear during the pulse, and levels 
outside the photoassociation window may remain pop- 
ulated after the pulse. 
— Another situation occurs when h{SLLj + IxItc) > 
since it is possible to transfer population to the con- 
tinuum, or to highly excited vibrational levels in the 
excited potential Ve (R) ■ The results of the calculations 
indeed show evidence for population transfer at large 
internuclear distances. After the pulse, due to the at- 
tractive character of Ve{R), with —C^/R^ asymptotic 
behaviour, the wavepacket created in the excited state 
is moving towards shorter distances R. For the up- 
permost levels, with an outer turning point located 
beyond ^ 500 ao, the vibrational half period T^n,/2 
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becomes comparable to the radiative lifetime, so that 
spontaneous emission may take place before the inter- 
mediate distance region is reached, the photoassoci- 
ated molecule decaying into a pair of atoms. In con- 
trast, the wavepacket corresponding to lower levels has 
time to reach the intermediate region where sponta- 
neous or induced emission may populate bound levels 
in the Vg potential. The present calculations indicate 
that besides the acceleration due to the potential, the 
laser pulse has given a "kick" to the wavepacket in the 
],Ve curve, thus reducing the time necessary to reach the 
intermediate distance region where radiative stabiliza- 
tion of the molecule may take place. Optimization of 
the pulses in view of increasing this "kick" is a promis- 
ing direction for future work, since at large distances 
the continuum wavefunctions representing the initial 
state display a large amplitude , and many vobrational 
levels of the Ve potential are located close to dissoci- 
ation limit, making the photoassociation process very 
efficient at small detunings, while it is well known from 
cw experiments that the bottleneck in such situation 
is the stabilization process ^ 

Both the present paper and Ref. [201 have been investi- 
gating the population transferred to the excited state due 
to photoassociation with one chirped laser pulse. In or- 
der to get final conclusions more useful to experiments, 
the theoretical work should develop further in two di- 
rections . First, besides the photoassociation step, future 
calculations should investigate the efficiency of the sta- 
bilization process, bringing molecules to bound levels of 
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the ground or lower triplet state, via spontaneous or in- Indeed, the classically allowed domain for a continuum 
duced emission. In particular the relevance of two-colour wavefunction is i?t < i? < i^, where Rt is the inner turn- 
experiments where a second pulse with a larger central ing point, and, in the standard semiclassical WKB form, 
frequency is transferring population to the desired levels the wavefunction can be written: 



should be analyzed. Second, although the present estima- 2. 

tion for photoassociation with a realistic repetition rate is yH^) ■> Rt 

promising, it does not take into account the fact that the with the Bohr-Sommerfeld quantization: 

second and further pulses are operating on a modified ini- /"^^ , / „/x ,„/ / 1 - 
tial continuum state, where population has been extracted 



^a.E=EAR) = iV„^==sin / {k{R')dB' + -) (53) 



/ " k{R')dR' = {n+]-)TT (54) 
JRt 2 



to be transferred to bound levels in the excited or ground The normahzation factor (iV„) is proportional to the 
state, or redistributed in the neighbouring continuum lev- density of states. 



els. The evolution of the atomic sample in the presence of , ,2 ^/^ 1 ^ /"^^ 

a sequence of short pulses is an important issue, particu- 
larly in view of possible applications to condensates. 



Therefore: 



^s,e=eAR) = \E=EX"^4>a,E„{R) (56) 

9 Appendix A. The energy normalization of 

the ground state continuum wavefunctions 10 Appendix B. Chirp rate in the time 

calculated in a box. domain for focussing the excited vibrational 

wavepacket at the inner turning point 

The s-wave ground state continuum wavefunctions (nor- 
malized in the energy scale: < ^q.eWq.E' >= S{E — E')) The gaussian linearly chirped laser pulse excites succe- 
have the following asymptotic behaviour: sively several vibrational levels in the molecular surface Vg, 

ro — • lhl'R\j- \ creating a wavepacket located at the outer turning points 

= V^^^^^^^^i^ (52) 

V Tin K{H) Qf corresponding vibrational wavefunctions. The cen- 

where r^g is a slowly varying phase and k{R) the local wave ter of the pulse excites sX t — tp the level vq. We choose 

number determined by the electronic ground state poten- the negative value of the chirp parameter % < in the time 



tial Vg{R): k{R) = l/hy/2ii{E — Vg{R)). For E — En, domain such as the excited wavepacket, located initially at 

the energy-normalized wavefunction ^g.E^ (R) is deduced the outer turning points in the 0^ potential, to be focussed 

from the unity-normalized wavefunction 4>g,En{R) using at the time t — tp + Tyib(vQ)/2 at the inner turning points 

the density of states ^ \e=e„ at the energy En jlH], 021 • of the excited vibrational functions. Then the chirp value 
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has to compensate the dispersion in the vibrational peri- 
ods Tyib{v) of the vibrational levels v resonantly excited 
at different times by the instaneous frequency u!(t)/2n of 
the laser. 

Let's note Ey the binding energy of the level v {Ey > 0, 
< 0). The levels vq + 1, vq and vq — I are reso- 
nantly excited at tp — (Ey„+i — Eyg)/{hx) < tp, tp, and 
tp — {Eyg-i — Ey^^)/{hx) > reaching their inner turn- 
ing points with the delays TyibivQ + l)/2, Tyn,{vQ) /2, and 
Tyib{vo — l)/2, respectively. We deduce the value of the 
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Fig. 12. The overlap between the initial stationary continuum 
chirp X from the condition that the components of the state in the a^r+potential, corresponding to Eo/kB = 54 fiK, 
wavepacket to be in phase at the inner turning point. In and the wavefunctions of the bound levels of the 0" potential, 



a first approximation the x is supplied by the condition: 



TyibiVo) 



E, 



V(, + l 



hx 



E. 



va-1 



TyibivQ + 1) 

2 

Tvib{vo - 1) 



nx 2 
Using the definition of the vibrational period: 
2TTh Anh 

Tyib(vo) = 



E., 



vo + 1 



(57) 



(58) 



with energies Ey between -6 and cm^^ under the 6s-|-6p3/2 
limit (the smaller value observed in the vicinity of -3 cm~^ is 
due to the tunneling from the outer well to the inner well of 
the 0^ potential). The horizontal arrow shows the "window" 
of Og vibrational levels between v—92 and v=106, excited by 
the pulse studied in Ref. |2U| . 



and of the revival period |44 |: 
47r?i 



^1 I 



£■110 + 1 + Ey„^l — 2Eyg 

one can deduce the following value of x from the Eqs. (|57|) 



corresponding to £'o/fcs = 54 fj,K, and the wavefunctions 
(59) of the bound levels of the 0^ potential, with energies Ey 
between -6 and cm^^ under the 6s-|-6p3/2 limit. 



X = -27r 



[Tvib{vo)]^ 



(60) 



11 Appendix C. Overlap of the initial 



continuum with the 0„ vibrational 



wavefunctions. 



Fig.ll2lrepresents the overlap integral | < 0„ Ey\'^S:^ Eq 
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